OTS PRICE 


s 

1 

A 



i N 6 j 1 9 9 4 5 

i 

A Contribution to the Theory of Critical Inclination 
' of Close Earth Satellites 

2. Case of Asymmetrical Potential^ 

| Shinko Aoki 

| Theoretical Division 

I Goddard Space Flight Center, NASA 

I Greenbelt, Maryland 

ABSTRACT 

| The previous theory (1962) is extended to the case of an 
asymmetrical potential of the main body especially in the 
case of small eccentricity by assuming that J n , the coeffi- 
cient of the n-th order zonal harmonics of the potential, is 

of n-th order of magnitude and the eccentricity is ofthe first 

order of magnitude. Here there is a peculiar k.nd of libra- 

tion which never occurs in the case of moderate eccentricity. 

This peculiar kind of lib rat ion splits into two kinds of 
libration, depending either on the anti-asymmetrical terms 
prevailing case or on the symmetrical terms prevailing case, 
which is a continuous transformation ofthe type described 
in the previous theory. Numerical test disclose s that for the 
earth the former peculiar kind of libration occurs. Also it 
is shown that the fifth coefficient is comparatively large and 
plays an important role in the asymmetrical theory for the 
case of smalleccentricity. The present paper shows that the 
anti-asymmetrical terms can not be neglected for the earth 
in the case of small eccentricity. 


tlhis work wos performed under o National Academy of Sciences Post-doctoral 
Resident Research Associateship program connected with NASA. 
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§ 1 Introduction 

In the previous paper by the same author (1962) the motion of a 
close satellite in the vicinity of the so -called critical inclination was studied. 
However, at that time the author omitted the effect of the anti-symmetrical 
terms in the potential o£ the main bodies. The theory Including these 
terms was given, for example, by Y . Kozai (1961). However, in spite 
of his comment on the case of small eccentricity, it does not seem 

very extensive. 

Therefore the author decided to develop a theory to avoid the 
difficulty connected with small eccentricity. A, preliminary consid- 
ration shows that this difficulty arises only in the vicinity of 

e = 0, and that there is no trajectory in general extending from the vicinity 

of e = 0 to a value of the order of the unity. Keeping this in mind, 
the author has expanded the Hamiltonian into a power series in 
e, which is assumed to be of the first order of magnitude, assuming that 
T the coefficient of zonal harmonics of n-th order, is of the n-th order 

J n * 

of magnitude, and that the deviation given by the following formula 
is assumed to be of the second order. 

a = 1 - 5H 2 Ora o )* 1 , (l- 1 ) 

where H, a constant, is the projected angular momentum to the equa- 
torial plane, and a o the mean semi-major axis. 

The necessary terms up to the sixth order of magnitude in this 
respect are picked up from the Hamiltonian F*. which is given by 
removing the so-called periodic perturbation terms. The terms 
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lower than the sixth order are all constant, therefore the sixth order 
terms are considered as the leading terms in our theory. 

The author found that the equation of motion can be reduced to the 
same type as the ones given by H. Andoyer: 

dh = dk _ (1 

dr 3h * dr 3k 


where 

T = - /3h + (7 + /3')h 2 + (y - /3')k 2 + (h 2 + k 2 ) 2 , (1.3) 

with a different restriction on the sign of /3' than that given by Andoyer. 
He assumed a restriction such as 

P > 0 , 

which is always confirmed in his theory connected with the librational 
problem of asteroids. In our theory, however, this restriction should be 
removed; namely, in some case we have P < 0 as well as p > 0. Espe- 
cially in the case of the earth, p is negative. Therefore, some change 
from his theory must take place. This slight change might easily be 
overlooked, but as will be seen later in this paper, some alternations 
are required afterwards, if we want to have real expressions for the 
solutions . 

The main purpose of this paper is to develop the theory with this 
difference in mind, and exclusively to give the case for ft < 0 . §2 pre- 

sents a preliminary discussion on obtaining the equations of motion in 
Andoyer T s form. §3 is the same as that of his theory; however, in 

order to avoid some confusion the present author rewrites the results 
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in a very compact form. In §4 some changes from his own appear in 
separating the several cases connected with the relation between the 
quantities y, = y//3 2/3 and /S' = P /P 2/3 . §5 gives a classification under 

which the solutions connected with the quantity u, = u//3 473 , where u is 
the energy constant of the system, should be written down separately. 

§6 gives the whole expressions of the solutions in the real number repre- 
sentation. §7 is an appendix which gives the characteristics, which 
means the trajectory within the plane of e cos g and e sin g without any 
attention to the relation with the time, the independent variable. Several 
pertinent numerical results are given. §8 gives some discussions related 
with the convergency, the relation to the previous theory and soforth. 
Finally §9 gives the conclusion, 

In order to test the assumption imposed on J n , the quantities for 
the earth are listed below - not only in actual values but also in units 
of proper powers of a 0 j3 l/3 = 0.0528, which is assumed to be of the 
first order of magnitude: 

J 2 = 1082.36 x 10‘ 6 = 0.3882 /S 2/3 a o 2 , 

J 3 = -2.566 x lO' 6 = -0.01 74 /3 3/3 a o 3 , 

J 4 = -2.14 x lO ' 6 = -0.275 /3 4/3 a o 4 , 

and 

J s = -0.063 x 10' 6 = -0.154 /3 S/3 a o s . 

Roughly speaking, these numerical values show that J 2 , J 4 and J s play 
approximately the same order of role in the vicinity of e = 0, say 
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e =0.10 or 0.05; on the other hand, J 3 does not take any important role 
here. 

2. Equation of Motion 

The potential, under which influence a negligible small mass parti- 
cle moves, is assumed in the following form: 


V 



1 



P n (sin 8) 


) i 


( 2 . 1 ) 


where /x = Mk 2 , k being Gaussian constant and M being the mass of the 
main body, 

r is distance of the particle from the center of the 
main body, 

Sis the declination of the particle, 

J n 's represent numerical constants which characterize the 
spheroidal potential; 

in this paper the summation 2 extends from 2 to 5 as seen later. 

Using a result of the so-called secular and long periodic parts in 
the original Hamiltonian such as given by Y. Kozai (1962a), we may 
pick up only the following necessary parts provided that J n is assumed 
to be of the n -th order of magnitude as well as the eccentricity e to 
be of the first order: 
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F' - — - — — (e 2 - a ) 2 

16 L 6 


+ 21 ^ J 4 e 2 


40 L 10 


^ c 2 cos 2g 

40 L 10 

3 J 4 2 9a 
e cos 2g 

40 L 10 


3 J 3 , . 

25 e(e 2 - a) s in g 


4 fS L 8 

3 M- 7 J 2 J3 
4~V§ L 


12 


e s in g 


M 7 J 5 


9 ^ J 5 

+ — — — e s in g, 


10 fS L 12 


( 2 . 2 ) 


where a = 1 - SH 2 ^)* 1 = aconstant which is assumed to be of the secondorder, 
and the other notations correspond to Delauney's. In this expression, the 
terms of beyond the sixth order of magnitude are neglected. As the Hamiltonian 

has only sixth order of magnitude and nothing else, neglecting the higher 
orders, we may take and 77 as the canonical variables, 


/Le sin g = Vl e cos g - V, 


in the canonical equations of motion: 

d<f _ (2.3) 

dt d Tj ’ dt 
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For mathematical simplicity, if the canonical variable s and the inde - 
pendent variable are changed: 

k = e cos g, h, = e sin g , 


r = 


3_ M 4 J 2 . 3 M 

16 L 7 


1/2 h 


t = — 


16 7/2 

a o 


t , 


(2.4) 


then the equations of motion will be 


dh i _ 3F i dk _ 3F i 
d t 3k ' dT 3hj 


(2.5) 


where by neglecting an unnecessary constant in F*, 

F = — — — (F* - constant) 

1 3 M 4 J 2 

= -^h + ( 7l + /S,')^ 2 + (7! - fli ) k 2 
+ J3 (h^ + k 2 ) hj + (hj 2 + k 2 ) 2 , 


with 


% 


4 J, _4_ £3 J24_ Js 

71 J 2 a o ‘ 7S a 0 3 " SYS J 2 a/’ 


\ 



2a, 



4 J 

/ 5 J 2 a 


> 


(2.7) 
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la order to remove the third degree terms in h, and k , h, is changed to h 
given by 


h = hj 



( 2 . 8 ) 


Further, if the equations are transformed in terms 


of p and q , 



(2.9) 


the equations can be written as follows. 


with 


dp 3F dq 0 - 

— = -2i — • t: = +2i — 

dr Bq dT a P 


( 2 . 10 ) 


f = - - /3(p + q) + rpq + \ fi 1 (p 2 - q2 ) + p2 q2, 
2 ^ 


( 2 . 11 ) 


whe r c 


and 


P - *7 + c 


2 Ja 2 

18 V 

16 ^3 J 4 

' 5 a J 2 a 2 + 

J 2 o 

25 J 2 a o 4 

25 J 2 3 a o 4 

24 J 3 Js 

3 J/ 

= cons t an t ^ 

25 J 2 2 B 0 4 + 

25 J 2 4 a o 4 

8 J4 2 J 3 2 

\ 4 J 3 

24 Js 

*j 2 2 s j 2 3 ( 

jT5a o 3" 

5/5 J 2 a o 3 ’ 

14 J 4 4 

J 3 2 



5 J 2 * 0 2 


5 j»- a 


2 o 2 


' 2 o 


/ 3 ' = 


5 a 


, j. Jji, 
J2 + j, ' j 2 2 1 


> ( 2 . 12 ) 
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The form of the Hamiltonian (2.11) shows that it is the quadratic 
in p or q, respectively; this was the technique which H. Andoyer (1903) 
used for the problem associated with the libration near the commensur- 
ability between the mean motion of asteroid and that of the disturbing 
body, Jupiter. Therefore, we can follow after his development in order to 
solve the equation of motion. Nevertheless, there is a slight difference 
between his case and ours concerning the sign of some coefficient. He 
assumed that the coefficient p is always positive; but in our case p 
can take a negative value as well as a positive one depending on the inter- 
relation of the magnitude of J n ’s. Especially for the earth it is negative 
as seen later. Therefore, we are restricted to take the negative P case 
here unless otherwise mentioned, because for the positive case any 
alternation from his results is not required. 


§3. General Expression of Solution 

The general form of the solutions does not change from Andoyer's 
results. But, in order to avoid some confusion about the real/complex 
criterion, we shall make a slight change in the notation. The calculations 

are omitted but only the results obtained are given here. 

The solutions of equations (2.10) are given in the following form: 


i Va 


i Va L 


- | Up + x) - 5 Up - X) * 

2 * J 


(3.1) 


C[t + X + -P + -> 


- — C(P + X) “ 2 ^ 
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with 

S<p) = -( 3l 2 +a 0 ' ajXao'V 1 (real), 

p '(p) = i(a' 2 a 3 + 3a 0 ' ai a 2 + 2 ai 3 )(a 0 ')' 3/2 (imaginary), 

p(X) =(P(P ) i ( a o ‘^ /2 P ' 1 W ( real ). 

p J (X) = -/S' /8‘ 1 (a 0 ') 1/a »>'(P) (imaginary), 

whe x' e 

a 0 '=-8/3\ a x = -2/3 , 

a 2 = i (8u + 2/?' 2 - 2y 2 ), a 3 = /3(y - P ), 
o 


a, =8u/3' - /3 2 ; 


u is the energy constant such that 


F + u = 0. 

jp(z) is the Weierstrassian elliptic function with the parameters, 

g 2 .= -a 0 ' a 4 - 4a t a 3 + 3a 2 , 

g 3 = 2SJ a 2 a 3 - a Q ' a 2 a 4 - a t 2 a 4 
+ a p a 3 2 “ a 2 ’ 


(3.2) 


(3.3) 


(3.4) 


(3.5) 
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such that 


p'(z) 2 = 4^ 3 (z) - g 2 tf>(z) - g 3 - 


v^.6) 


where p'(z) denotes the derivative respective to z. C(z) is the associated 
£ -function defined by 




(3.7) 


In the expression of solution (3.1) X is an integration constant, which 
with the energy constant u forms a system of the arbitrary constants in 
the solution of the equations of motion, u is of course a real value; X is, 
however, not necessarily real. X should be taken such that the expression 

(3.1) could give real values of h and 4 k. 


h = 1 (p + q), 

k = j7 (p " q >‘ 


(2.9)’ 


§4 


Discrimination Among Several Cases 

At the first step it is necessary to know the sign of the discriminant 


A, 


A = g 2 3 - 27 g 3 2 . 


For brevity, let 


p = p(x) = - (4u - 7 2 + 6/3*7 - 5/3' 2 ), 

3 

and 

Q = p'( X )/i ^ 16u/3' - 4/3' (/3* - 7) 2 - 2/3 2 , 


(4.1) 


"N 

>- (4-2) 

J 
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which are also real values even for our case, then 


g, = 12P 2 -4(y- 2/3' ) Q , 


and 


S3 


= -8P 3 + Q 2 + 4(7- 2/3' )PQ. 




It follows accordingly that 


A = 1024Q 2 J u 3 - i (7 + fi' ) 2 u 2 


+ _L (7 + /3' ) [(7 + /3 ') 3 + 9^ 2 ] u 

16 


1 Q2 


64 


/3 2 


(7 + /3') 3 + 


For brevity, let 


u, = u//3 


4/3 


/3.' = /3'//3 2/3 , 


7. =7//3 2/3 , A, =A// 3 8 , 


Q. = Q//3 2 , P. = P//3 473 , 


and so forth, then 


A. 1024Q, 2 ju . 3 -i ( 7 . +/3.')u , 2 + ^ ( 7 . [( 7 . + P'.) 3 + 9j 


64 


(7. + /5;) 3 + X 


(4.3) 


(4.1)' 


(4.4) 


u . 

(4.1)" 


If the subscript * is omitted in this expression, then every term remains 
the same as in (4.1)' except for fi which is replaced by the unity; besides, 
as easily seen, fi is of the third order of magnitude, and we may 

choose 


h. = h // 3 l/3 , k, = k/ fi l/3 ■ 


(4.4)’ 


As a result, all quantities are measured in units of proper powe 
of which is assumed to be of the first order of magnitude, provided 

that d 7 0. Therefore without any restriction we may take d 
unless otherwise stated. This is an application of the non-dimensional 

analysis . 

in any way, there are two cases for A = 0. 


i) 


Q = 0- 


(4.5) 


If we consider (4.5) as an equation for u , then we know the value u Q which 
satisfies this equation*. 


U„ I ' r ' 1 - 


(4.6) 


ii) 


f(u) 


3 1 / p, ‘ \ , , 2 4. Z-±JL- [(y + fi') 3 +9]u 

u) s u 3 - x (7 + fi ) u + 16 L, ' / ^ ' 


l . 3 22 _ o. 

- 64 + P ) - 4 ' 


(4.7) 


divided according to the sign of the discnmi- 


Here two cases again are 
nant A # of f (u) itself. 

=-i5s {(’••’’At) 1 


(4.8) 
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namely, if 


y + fi' > -§. 


(4.9) 


then among the three roots Uj, u 2 , u 3 of Eq. (4.7) in terms of u, two of them must 
complex values, while the other, say u x , remains real. On the other hand, 
if 


7 + /S' < 



(4.10) 


then the three roots are all real, 

U 3 <U 2 <U 1. 

Now, it is necessary to know in what sequence the roots u. (i = 0, 1, 
2, 3) appear; in other words, to know the relation between u Q and 
u. (i = 1, 2, 3) (A ' > 0) or u x ( A' < 0). In the case when (4.9) occurs we 
can easily determine it in the following way: first of all we have 


f(u) < 0, 


i f and onl y i f u 


> 

< 


u , 


however, since 


f(V=£(r ! 


♦ 5 * 


I* 1 


>P f 


(A) 

(4.11) 


where P Q denotes the value of P when u = u o , we have, if fi' > 0, the second 
hand side of (4.11) is always positive, but if ft' < 0, this is not always true. 
Consequently , we have a new discrimination, which never occurs in the 
case when /?'> 0: 

if y 2 +2 z 8 '” 1 > 0 then u o> u i‘ ( B ) 

On the other hand, in the case when (4.10) is satisfied, we may only say 
that 
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if 7 2 + | /3" 1 < 0, then u 2 <u 0 <u,, 


or 


U 0 < U 3 ; 


(C) 


and if y 2 + | /3'" 1 > 0, then 


u i X u o' 


or u 3 '< u Q < u 2 


The detailed criterion will be given in a later part. 

In the second step it is necessary to divide the cases according 

to the sign of (g 3 ). (i = 0, 1, 2, 3), where (g 3 ) i denotes the value of g 3 

when u - u . This is important because, for example, if (g 3 )i < 0 then 
i 

the energy constant u. corresponds to the unstable equilibrium point (s). 
It is easily confirmed for (g 3 ) 0 : since 


(g 3 >o = “ 8P o 3 




(4.12) 


when J3 1 < 0 , 

if y - + | /3'~ 2 > 0- then (g 3 >0 > °- (°) 

In order to know the relation between u Q and u. ( 1 = 1, 2, 3) in more 
detail as well as the sign of (g 3 ) it it is convenient to divide the y vs. /?'- 
plane into several parts so that in each of them the situation never changes. 
For this purpose, it is sufficient to draw lines where the situation would 
change. The following four lines serve this purpose: 

y + /3'=-|. (4.13) 
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y + (-2 /5' )- 1/2 - 0 

(4.14) 


7 _ (-2 ,6')- 1/2 = 0 

(4.14) 

and 

7 +1 = 0. 

(4.15) 


Iv. effect. u Q = u. occurs only on the lines (4.14) or (4.14)', u 1 - u 2 or 
u 2 ,u 3 only on the line (4.13); and (g 3 ) 0 = 0 only on (4.15), (g 3 )i= 0 

(i 0) occurs on either (4.14), (4.14)' or (4.13), as is easily confirmed. 

Before dividing the plane by the four lines given above, for the sake 
o" simplicity let us consider the division by separating the negative fi 
h df-plane into three parts: 

Case!. ^'<-2, 

in this case we have 

_ |> + (-2/3')- 1/2 > - ( -2/3') _1/2 > /3' - | (,/3')- 2 . (4.16) 

Case II. -2 < /3 < - 2 > 

(-2 / 3')- 1/2 > - Z 3 ' > -(- 2 / 3') -1/2 > F ( Z 3 ') -2 - < 4 ‘ 17 > 

Case III. ~ \ < Z 3 ' < °> 

( - 2 /Q') -172 > -/S' > - (-2/3')- 1/2 > fi' -\U$')- 2 ■ (4 ' l8) 
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The reason why we make a discrimination between Case II and Case III 
is that 

h then 

-/S' -! = -(-2/3')- 1/2 =/r-|/3'- 2 = -l. 

so that a domain, -/?' - | > 7 > - (- 2 ^')‘ 1/2 for example, is cut 
out at /3' = _i and there is no continuous route connecting the regions 
corresponding to Case II and Case III without meeting any one of the 
above four lines. 

Now, we sub-divide the regions with respect to the four lines 
stated above; for example, Case I, means fi' < -2 and 7 > " fi - tj, Case 
I and - / 3 ‘ - | > 7 (-2 /3 ' ) ~ 1 7 2 and so forth. This sub-division 

is not essential 1 as a whole because Case I, , II, and HI,, for example, 
have no differences in the sense of continuous deformation. Therefore 
it is useful to assembly some of the too sub-divided cases into one 
case, then finally we have the following division: 


Case 

i 

which involves 

cases 

I 1 * n i ’ m i 5 

Case 

2’ 

which involves 

cases 

V 

Case 

2" 

which involves 

case s 

n 2 , m 2 ; 

Case 

3’ 

which involves 

cases 

C n 3 ; 

Case 

3" 

which involves 

cases 

HI 3’ 

Case 

4‘ 

which involves 

cases 

l 4 ’ U 4 

Case 

4" 

which involves 

cases 

nI 4 ; 

Case 

5 

which involves 

cases 



1 There is, nevertheless, a practical advantage, because /3', or speaking more precisely 
J3> # = fi' // 5 2/3 , depends only on the coefficients of the zonal harmonics for the potential 
and is totally independent of the initial condition. 
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The Fig. 1 shows the division in the /3' vs 7 -plane. Here for the 
sake of comparison, the case where (5' > 0 (Andoyer’s sub-cases <a , 
a* , b', and b" ) is added. 

We evaluate u, and (g.). in each case either by giving extreme 
values for /3‘ and y or by giving special numerical values as 
given’ in section 7. In any case one set of valued for /3' and 7 
is enough to determine the sign of (g 3 ). an d relation of u 0 
and u . ( i = 1 » 3). 

As an example, the case 3' is shown where extreme values are given. 


Case 3 1 


/ 3 1 = - £~ 1/3 , y = + , e — + 0 , 


and 


u - i e~ 2/3 + - e 1/3 

u o - 4 8 


/ A 512 a 

^3) 0 - + 27 e 


u - - e - 2/3 ± — e- 1/6 , 
u i,2 - 4 e 2 


„ . _ 512 _ 2 

(£ 3 ) 1.2 " + 27 e * 


U.=-J «‘ /3 . 


, . _ j_ - 2 

(£ 3)3 28 6 


where unnecessary higher order terms are omitted. 
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Therefore 


and 

(g 3 ) 0 > o, (g 3 ) i > 0, (g 3 ) 2 > 0, and (g 3 ) 3 < 0. 

Table 1 shows the results. From this table, combined 
with (4. 7), it is easy to find the sign of A(u) (e .g. , in the case 3 1 , if u 1 <u 
then A(u) >0,ifu Q <u<u 1 oru 2 <u<u Q then A(u) < 0 , and so forth). 


Table 1 

Roots of A(u) = 0 and Sign ofg 3 at Respective Root s 


r r 

j Cases 

U i 

( g 3 ^0 

Cs 3 >i 

(g 3 ) 2 

(g 3 ) 3 

I Ca s e 1 

I 

U 1 < u 0 

r 

T 



i Case 2 ' 

\ 

U 3 < U 2 < U 1 < U 0 

+ 

+ 

+ 


| Case 2 * 

u o < u i 

r 

+ 



| Case 3' 

1 U 3 < U 2 < U 0 < U 1 

+ 

T 

+ 

- 

| Case 3" j 

U 0 < U 3 < U 2 < U 1 

+ 

+ 

- 

+ 

Ca s e 4' 

U 3 <U 0 <U 2 <U 1 

+ 

+ 

T 

- 

Case 4" 

U 3 <U , 0 <U 2 < U 1 

+ 

4. 

- 

+ 

Case 5 

U 3 < U 0 < U 2 < U 1 

- 

+ 

+ 

+ 


§5 Values of p and y 

For the purpose of giving a real expression of the solutions, it is also 
necessary to have the range in which p ory falls. As a first step, a com- 
parison is made here between #>(p) or ^(x) and the parameters of ^-function, 
namely e x , e 2 , and e 3 . Let 

4 P 3 (z) - g 2 P ( z ) - g 3 = 4 O(z) - e !> (s»(z) - e 2 ) (p(z) - e 3 ) , (5.1) 
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then it is well known that if the discriminant A > 0, the three parameters 

e , , e , , and e , are all real; contrary to this , if A < 0, then only one of 
12 

them, say e x , is real. 

Since 

p'O ) 2 =4 (p(p) -ej) (fp(p) -e 2 ) (p(P) -e 3 ), (5>2 

P(P)=-C aj 2 +a 1 2 a 2 ) a 0 ' is real, ( 3 - 2 


and 


'(p) = i (a 0 2 a 3 + 3 a Q ' a 1 a, + 2a 3 ) (a 0 ')" 3/2 is imaginar 


S 3 


it follows that, when A > 0: 


p(p) < e 3 or e 2 <p(p)<e l , 


(5.3) 


where the three arguments are denoted such that 


e 3 < e 2 < e i 


Similarly, with 


£>(X) = f?(P) +4 


y P" 1 tp'(p) - p is real 


1 




(X) =-/!' p ' (p) s iQ is imaginary, 


(4.2) 


when A > 0: 

p(X ) <e 3 or e 2 < s?(X) < e r 


(5.4) 
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In the case when A < 0 , it follows necessarily that 


*»(P) < e i 


and. 


(5.5) 


s>(X) < e i • 




Besides, from (4.2) it follows that 


>(P) > p(x)< Q > 0 


that is to say (/?' < 0), 


u < u. 


(E) 


The above is simply obtained. The next is to decide whether *><X) and 
p(p) drop either in the region between e x and e 2 or the region smaller 
than e 3 . This decision can be given by knowing the behavior near u = u Q , 
where 'generally the situation is changed, and by knowing the behavior in 
the extreme case such as u - + oc. The detailed calculations are omitted 
here, and only the results thus obtained are listed in Table 2. In this 

table, for brevity, we let: 


p = fp(x) and (p = s?(P) • 
Besides, <p '(/ )/i andja ' (p)/i are determined by 


(5.6) 


fP 


(X)/i = Q = 16/3'(u-u 0 ), 


(5.7) 


and 


& 


(p)/i = - 16(a')~ 1/2 (u - u 0 ) 


(5.8) 
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which give us the sign of both values, provided that (a 0 ') 1/2 » or more 
rigorously (a 0 ')* /2 s (a ' 0 ) 172 //3 173 , is non-negative without any restric- 
tion. In effect, since in the form of the solutions h and k , or more 
rigorously h* = h//3 1/3 andk* = k//3 173 have (ap)“ 1/2 , or more rigorously 
(ag)iT 172 = /3 173 /2(-2/3') 172 , as a factor, the change of sign of (a^) 172 
as a whole does not produce any real difference in the solutions at all. 
It is to be noted that when /3<0 this provision means (a^) -172 itself, 
which does not mean a non-dimensional quantity (ag)£ 172 , should be 
taken as negative. By this method we may ignore the difference 
of sign of j3 but may unify the procedure as a whole. This is justified 
directly from the fact that in the equations of motion the sign of /3 can 
be changed without any significant alternation except for the change 
of signs of h and k for the same time. 
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Table 2 


The Values of P = &»(X) an< ^ P ~ p(P) (P < 


Cases 

Range of u 

Values of P and p 

Sign of 

and^' (P) 

i i 

Classi f ication 2) 

Case 1 

u 0 < U 

P < e 3 < e 2 < P < e i 

- 

01 


u = u o 

P = e 3 = e 2 = P 

0 



Uj c u < u 0 

p < e 3 < e 2 < < ej 

+ 

02 

i 

U < Uj 

P < P < e, 

+ 

12 

Case 2' 

u 0 < u 

9 < e 3 < e 2 < P < e i 

- 

01 

1 

J 

i 

U = U 0 

P ~ e 3 ~ e 2 ~ P 

0 


1 

1 

U 1 < u < u 0 

! P < e 3 < e 2 < ® < e i 

+ 

02 

j 

j 

u 2 < U < Uj 

P < P < ej 

+ 

12 

i 

| 

u 3 C u < U 2 

P < P < e 3 < e 2 < e i 3) 

+ 

• 22 

i 

1 

| 

u < u 3 

P < P < e ! 3) 

+ 

12 

; Case 2" 

Uj < u 

fP < e 3 < e 2 <P < e x 

- 

01 

! i 

U 0 < u < u ! 

P < P < e i 

- 

11 

} 

u = u 0 

P = e 3 = e 2 = P < G 1 

0 


| 

4 

u < u 0 

P < p < e i 

+ 

12 


2) See Table 3. 

3) e x in the lower line is the continuation of e 3 in the upper line as a result of the dropping of e t 
and e 2 in the upper line into imaginary values. 
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Case 3* 


Case 3 " 


Table 2 (Continued) 

Sign of 


Cases Rangeofu Values of P and p P^X) and Classification 


Uj < u 


P < e 3 < e 2 < P < e, 


U < u < Uj ip < P < ej 


= e„ =P 


u 2 < U < u 0 p < p < e i 

U 3 < u < u 2 p < p< e 3 < e 2 < e x : 

u < u, P <P < ej 3) 


Uj < u 
U 2 < U < U 1 


< e 3 < e 2 < P < Cj 


P ' c 3 2 

P<P< ej 41 


U 3 < U < u 2 p < P < e 3 < e 2 < e i 

U 0 < U < U 3 S 3 < P < e i 

ii — n r G-, — P ^ e. 


Case 4' 


u < u r 


u x < u 


P < P < e i 

< e < e 2 < P < e ! 


u < U < Uj I p < P < e i 


< u < u~ e, < e 3 


< P < e, 


u = u P = e 3 = e 2 -P 

2 2 

u, < u < u 0 P < p < e 3 < e 2 < e i + 

3 31 + 12 
u < u 3 P < P < e x J) J ^ 

_ 4) . ,h. upper line i. .he cen.inue.ion el e 3 in *. oe • ™»l* °< *• d ^ PP ‘" 9 ** 

and e 2 in the lower line into imaginary values. 


24 



Table 2 (Continued) 


Cases 

Range of u 

Va lues of P and 

Sign of 

«-'<X> and^lM 

1 1 

n 

Classification 

Case 4" 

u « < u 

l 

IP < e 3 < e 2 < P < e i 


01 


u 2 < U < Uj 

*><?< e/ 0 


11 


U 0 < U < U 2 

V < P < e 3 < e 2 < e x 4) 

- 

21 


u = u o 

P ~ ? “ e 3 ” C 2 

0 

j 

j 


u 3 < u < u o 

e 3 < e 2 < P < p < e ! 

-f- 

24 

I 


u < u 3 

P < V < e, 

+ 

12 

Case 5 

Uj < u 

IP < e 3 < e 2 < P < G I 

- 

01 

i 

| 

u 2 < u < U 1 

. 

IP < P < e x 

- 

11 

j 

u 0 < U < U 2 

e 3 < e 2 < S 3 < P < e i 

— 

23 

[ 

. 

u = u 0 

OJ 

II 

II 

PL. 

11 

<N 

0) 

0 

j 


u 3 C u C U 0 

e 3 < e 2 < P < jp < ej 

+ 

I 24 


U < u 3 

L — 

P < f> < ej 

+ 

12 
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The next step is to determine the range in which p or y will 
fall. For this purpose it is convenient to classify the various cases 
which are associated also with the value of u . This was done in 
the last column of the Table 2. The specification of each class is 
given in Table 3. Thus we can determine the range of p and x> ** 
we remember that: 


whe n A < 0 , 


if 0 < v < oo " , then 
-oo<£>(iv) < e and #> # (iv)/i < 0, 
but if cjj < v < 2 , then 

e 1 >£>(i2' / )>“°o> and £p l (iv)/i>0, 


\ 


> 




(F) 


where i is the purely imaginary semi-period given by 


CO 


u 

1 



dp 

- P) (e 2 - p) (e 3 - p) 


when e x is real. 

On the other hand, when A > 0 , 

if 0 < v < oj" t then 
e > p (a; f + i v') > e 2 , * (w # + i v) / i > 0 , 

- oo < fj(i v) < e 3 , and gp 1 (i v) / i <0, 
but if to" < v < 2cc a , then 
e 2 < £> (a> ' + i v) < e r ^ ' + i v) / i <0 

e 3 > #? ( i i') > - oo and p # ( i y ) / i > 0 » 


N 


> 


(5.9) 


(G) 
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Note: Class 0 and Class 2 are corresponding to A > 0. On the other hand, 
Class 1 is corresponding to A < 0. 



where co' is the real semi-period and «*i is the purely imaginary semi- 
period given by (e 3 < e 2 < ej) 


CO 




dp 

) (P - e 2 ) (P - e 



(5.10) 


and oj h 



dp 

y 4( ei -P) (e 2 -P) (e 3 -P) 


(5.11) 


§6 Real Expression of Solutions 

As seen in the preceding section the values p and X are uniquely 
determined within the parallelogram, which is (2 a' , 2w"i)for A > 0 
or (2 , co[ + i) for A< 0 ; however, the parameters which are in- 
cluded in the expression of solutions are not p and X respectively but 

1 p and ix. These values can not be uniquely determined in the parallel- 

2 2 

ogram but correspond to the four values respectively. In effect, v 
is uniquely determined byJ p(v) = given constant , p (i/) has given sign, 


then 


1 

-v + o; 


1 

-V +0) 


u 


1 , 


and + o>* 


+ co i 


for A > 0 , 


or 


V ^ CO. 


1 

~2 V 


1 , / ti • \ 

+ o-("l + W 1 O. 


. 1 3 / 

and y v + y"i 


1 „ 


i, for A < 0, 


are all required values as well as 2 • 
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This phenomenon, however, does not provide any essential difficulty. 

In effect, the change from-t v to ^ + a.'' produces only a half period 
advance, and this is cancelled out by choosing the proper value for the 
additive constant On the other hand, the change from|y to 
l v + «"i or to 5 ^ + 5 + W ») produces another series of solu- 

tions, which satisfy the equations of motion as well as the given energy 
constant u. Speaking from the results advancely. Class 2 (which includes 
Classes 21-24) has two series of real solutions; contrary to this, Class 1 
has only one series of real solutions while the other corresponds to 
imaginary solutions, and Class 0 has no series of real solutions. This 
change can be also carried, out by changing X. to A. + oj i 
y + 1 . + <x> " i ) without any alternation of and-j y respectively, 

as was done by Andoyer, since the solutions are periodic qua function 
of or-ty respectively, and the argument-is in the form of 

x +± p ±± - x . Accordingly, we may fix the values for and^-y- And from 


and 


j(p +X) =x lt 


> 


~2 (P - X) = X 2 


( 6 . 1 ) 


we may construct Table 4. 

Thus we have obtained the real expressions of the solutions 
in each class. From (3.1) the expressions of h and k are easily obtained: 
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Table 4 


Constants Xj and X 2 
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h =~ (p + q) 


= U(r + X + Xj ) - l(r + X -Xj) 

4 i V— 2/3 ' 

+ C(T + X +x 2 ) - C(T + K -x 2 ) 
- 4(2 Xj) - 4(2 X 2 ) ] 


( 6 . 2 ) 


k = jj- (p - q) 

= I [ 4(r + X +Xj) + 4 (t + x -X, ) 

4 y-2 /3 ' 

- 4 ( T + x +x 2 ) ~ £ ( T + x - x 2 ) ] . 

When A > 0, the expression of 4(v) is given by 


77 i 

2 a/ 



E <T .077 .... 

sin-^v + a; O. 
n= i 


(6.3) 


where the double signs should be taken according to the sign of the 
imaginary part of v : 
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upper sign for 0 < Im(v) < 2oj 


or lower sign for -2 < Im(^) < 0, 
q = exp [-t Ta)"/aS] , 

= £<«'). (6 - 4) 

V 1 is the real semi-period, and to" i is the purely imaginary semi-period; 
when A < 0 , 


S(v) = 


V,v _ 


rr 


+ — ? i 

2 


277 X ' ^1 n77 / _ 1 ' 1 // \ 

) sm— 7 v w t w l , 

L l-qf W 1 \ 2 1 2 1 / 


(6.5) 


where the double sign should be taken according to the sign of the 
imaginary part of v : 

a 

upper sign for 0 < Im(v) < oj ^ 
or lower sign for - < Im(v) <0, 


q, = exp 


r * O' 

to + co 1 1 


2 to* 


- 77 1 


( = iq 2 , q 2 being a real pos i t ive value ) , 


( 6 . 6 ) 


1 ?! = £(<^). 


The real semi-period to given by 


32 


, f dp (6.7) 

"l = ■ — * 

J V4(p - ej) (p - e 2 ) (p - e 3 ) 

and oj" i is the purely imaginary semi-period. (Note: to* and — +"i 0 

form a system of the fundamental semi -period.) 

All the necessarypreparations have been completedand we are now ready 
to construct the real expressions of the solutions in the following: (P < 0) 


A < 0 (Class 1) 


N 


h = 


j—jl [y q ^._ ic 

CiW L^-* 1 + ^ n 


COs5 (T+\) +c; 

W 1 


V" _lL_ jc; 

La 1 - q2" l 


cos (r + A.) + C 

“i 


"1 

no r 


+ A 


\ (6.8) 


k = 


_1 TT_ y _ 

-o 7P 1 L* x 


9 2 „ , . nw . x 

2 - S' sin— x (r + X) 






n o>: 


y -1 — — ? S " s in -^7- (t + A) 

La l_ q 2n n "l J 


/ 


where 


q = exp 


2 w' 


_ _jL ,,1/2 T a -7/2 t 
t = 16 m J 2 a o 
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S' means the summation extending only over positive odd numbers, and 
S" only over positive even numbers. The coefficients are expressed 


by 


n7r X i * " W X 2 

C ' = c, ' cosh — + cosh — 

n l jm* e, >1 




a) X 




HTT/j* . n7TX 2 
C" = sinh r + sinh — —r 

n t r At) - 


n77 X] 

S ' = s, sinh — 

n 1 a )' 


- sinh 


n77X 2 


n77Xj n7r X 2 

S " = s inh ; sinh 


w, 1 


> 


(6.9) 


/ _ O 

no 2 


C ' = —— cosh 


n7T Xi 1 u n7T X 2 

cosh- 


>i 2 




1 n77 Xi" < 1 . u nw x 2 

c " = - — s i nh : -Sinh- 


oo{ 2 






The coefficients and the arguments involved in the above expressions are 

given for each class in Table 5 respectively. 

A > 0 (Class 0 has no real solution; on the other hand Class 2 has 
two series of real solutions; the existence of two sub-classes corresponds 
to this phenomenon. The second sub-classes are constructed by making 
X = V + a/'i in the first sub-classes respectively. After being expanded 
into series, the unnecessary prime on \ is omitted in the following.) 
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_J__.IL \y -^(c„cos ^(T + X).C no } , 
Z_J 1 - q I *> J 

n “ 1 


+ A 


( 6 . 10 ) 


- _J_ _£ Y - /s B sin2?(r + A)}. 

- yT 2 ? «' Z_i 1 - q 2n l " J 

n ~ 1 


where 


77 CO 


exp 


16 


~M l/2 J 2 a o" 7/2 


The coefficients are expressed by 


nw/j" . ^ n7r X 2 

C = c, sinh r- +sinh — 

n 1 a/ W 


n7T Xi" , n7T X 2 

S = s. cosh ; cosh 

n a <v> 




( 6 . 11 ) 


1 nTT'/J" i . , n7T X 2 

C = - — sinh ; — - - s inh — — 


The coefficients and arguments involved in the above expressions are 
given in Table 6. These constants are chosen such that all the expansions 
should be convergent, or in other words, the restriction imposed on (6.3) is 
always satisfied. If this would not be satisfied, before being expanded into 
series the argument has been changed keeping in mind the periodicity of 
h and k qua function of Xj • 

It seems necessary to add the dimensionality of to, X x * X 2 > etc - They 
are all expressed in the following way: u, = a;/3 2/3 . and so forth. 
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§7 Characteristics 


Some numerical calculations connected with the characteristi 
are given here. The characteristics means the orbit without any 
regard to the independent variable r but only the plotted line in h 
and k space. In this problem it is easy to obtain such a characteristic. 

In effect, 

~ « (7 

F + u = 0 ' 

gives one; or, in more detail 

(h 2 + k 2 ) 2 + (7 + /S')h 2 + (7 - /3')k 2 - /3h + u = 0. (? 

The same procedure of the non-dimensional analysis taken in §5 gives 
the following form (omitting the subscript v ). 

(h 2 +k 2 ) 2 + (7 + /S')h 2 + (7-/3')k 2 - h + u = 0. (7> 


from which we have simply 


k 2 = - jh 2 (7-/?')} 


j/-2/3* h 2 + h - u + i (7 - P ? 


Accordingly, by putting 


U, = -2/3' h 2 + h + i (7 - /5') 2 


and 


U 2 = -h 4 - (7 - /3')h 2 + h. 


the following criterion is given; 


(7.3) 


(7.4) 
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i) if u > Uj , then 0 solut ion for k; 


11 


)j if u < Uj, h 2 + i (7 -/?')> 0 , 


ii) n and u > U 2 , then 0 solution, 


ii ) 12 or u < U 2 , then 2 solutions; 


ii ) 2 if u < U r h 2 + ]■ (7 - /S') < 0, 


ii ) 21 and u > U 2 , then 4 solut ions, 

ii ) 22 or u < U 2 , then 2 solutions. 


> 


(H) 


The critical values for u also come from these criterion as follows: 


u o = u i ( h o>- 


3h 


0 , 


h = h 


u 

1 


= U 2 (h.), 



= 0 (i = 1, 2, 3). 


(7.5) 


(7.6) 


As an example, one set of P>‘ and 7 for each (I ,...111 ) is 
taken and several characteristics corresponding to some values 
of u are given in Figs #i 2-16. Table 7 presents the numerical 
values of u. (i = 0, . . . 3) etc. Also added is the Case E where /?*' = -0.13, 
which is, of course, included in Case III, but this is the actual case for the 
earth. Adopted values of the harmonic coefficients are fromKozai (1962b). 
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> 


(7.7) 


J 2 = 1082.36 x 10~ 6 J 3 = -2. 566 x 10' 6 

J 4 = -2.14 =< lO" 6 J 5 = -0.063 x 10 -6 . 


Accordingly , 


(5 = 0.147 x lO -3 • a" 3 , 


/3 1/3 = 0.0528 a" 1 


/3‘ = -0.360 x 10" 3 • a q 2 
y - -2a. - 5.540 x 10 -3 • a" 2 


(a = 1 - 5H 2 /L 2 ) 


> 


y3" = -4. 24 x 10 -3 a" 1 ; 


therefore , 


= fi' /B 2/ 3 = - 0 . 129 
y, = - 71.7 aj a - 1.99 
h, = 18.9 a Q h =18.9 a Q e sing - 0.0201 
k, = 18 . 9 a 0 k = 18.9 a 0 e cos g, 





(7.8) 


(7.9) 


where a 0 stands for the mean semi-major axis, e 
and g for the argument of perigee. 


for the eccentricity, 


It 

mainly 


is noted that in the case of the earth the asymmetrical part comes 
from J s , but not from J 3 . The linear shift term (/?" ) which de- 


pends only on J 3 in our theory is small. It is also noted that 
dependent on the mean semi-major axis a Q but only on J n ’s; 
if the potential is given, it is an absolute constant. 


/?' is not 
therefore , 


40 


In concluding this section, we should add some remaiks 
on the critical value u Q . As is seen from Table 7, the critical value 
Uq corresponds neither to stable point nor to a unstable one except foi 
Case 5, where of course it corresponds to the unstable points. The 
sign of (g ) 0 is positive except for Case 5. Therefore it might be con- 
sidered that u Q corresponds to a stable point. In fact, however, it cor- 
responds to imaginary equilibrium points. Or in other words, if we 
draw the characteristics corresponding tou Q it is split into two circles 
given by 



and for Case 5 two circles meet at two points so that the two points become 
unstable points; on the other hand, except for Case 5, they cannot meet 
each other or the radii(us) of one or both circle (s) become (s) imaginary. 
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Table 7 

and Corresponding h. and k 


Critical Values of u* 


Case I /3' = -3.0 


Sub- case 

7. 

i 

U *i ! 

h *i 

k *i 

. Remark 

i. 

+ 2.0 

0 

+ 6.2083 

“ 

- 


l 


1 

+ 1.0548 

+0. 8847 

0.0 

s t ab 1 e 

i. 

+ 1.0 

0 

+ 3. 9583 

- 

- 




1 

+ 2.0562 

+ 1.1072 

0.0 

stable 



2 

+ 0.0734 

-0.8376 

0.0 

stable 



3 

_ 0.1295 

-0.2696 

o 

o 

uns t ab 1 e 

i 3 

0.0 

1 

+ 3.5139 

+ 1.3008 

0.0 

s t ab 1 e 



0 

+ 2.2083 






2 

+ 1.0702 

-1.1309 

0.0 

s tab 1 e 



3 

- 0.0840 

-0.1699 

0.0 

uns t ab 1 e 

i 4 

-2.0 

1 

+ 7.8553 

+ 1.6290 

0.0 

s t ab 1 e 



2 

+ 4.6947 

-1.5286 

0.0 

stable 



0 

+ 0.2083 

| - 

— 




3 

- 0.0500 

-0. 1004 

0.0 

unstable 

i s 

-4.0 

1 

+ 14. 1383 

+ 1.9056 

0.0 

stable 



2 

i 4-10.397 6 

-1.8340 

0.0 

s t ab 1 e 



0 

+ 0.2083 

-0.0833 

±0.7022 

t 

uns t ab 1 e 



3 

_ 0.0358 

-0.0715 

O 

d 

stable 


42 



Table 7 (Continued) 


Case II J3, = -1.0 




■ 

u *i 


m 


I 1 ! 

+ 2.0 


+2.1250 



— 




1 

+0.2148 

+0.3855 

0.0 

s t abl e 

“a 

0.0 

1 

+ 1.0548 

+0.8847 

0.0 

stable 



0 

+0.1250 

— 

- 


IJ 3 

-0.6 

1 

+ 1.6022 

+ 1.0221 

0.0 

stable 



0 

-0.0850 

- 





2 

-0.1524 

-0. 6396 

0.0 

stable 



3 

-0. 1698 

-0.3823 

0.0 

uns table 

1 1 4 

-1.0 

1 

+ 2.0562 

+ 1.1072 

0.0 

stable 



2 

+0.0734 

-0.8376 

0.0 

stable 



0 

-0.1250 

— 

- 




3 

-0.1295 

-0.2696 

0.0 

un s t ab 1 e 

ri s 

-2. 0 

1 

+3. 5139 

+ 1.3008 

0.0 

s t ab 1 e 

l 



2 

+ 1. 0702 

-1.1309 

0.0 

s t abl e 



0 

+0.1250 

-0.2500 

±0.6614 

uns table 



3 

-0.0840 

-0.1699 

0.0 

stable 

i . 
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Table 7 (Continued) 


Case III ft', - -0.25 


Sub-case 

HI 

i 

u .i 

h *i 

— 


m. 

+ 3.0 j 

0 

+ 0.7656 

- 





1 

+ 0.1371 

+ 0.2645 

0.0 

stable 

m 2 

0.0 

1 

+0.5825 

+ 0.6959 

0.0 

stable 



0 

-0.4844 




hi. 

-1.3 

1 

+ 1.5506 

+ 1.0111 

0.0 

stable 

3 


1 

2 ; 

-0.1717 

-0. 5966 

0.0 

uns t ab 1 e 

1 

\ 


3 1 

-0.1777 

-0.4145 

0.0 

stable 



0 

-0.2244 

— 

" 


III, 

-1.5 

1 

, +1.7641 

+ 1.0546 

0.0 

s t ab 1 e 

4 


2 

| -0.0815 

-0.7296 

0.0 

uns tabl e 



0 

| -0.1094 

| 

— 




3 

-0.1512 

-0.3249 

0.0 

stable 

III. 

-3.0 

1 

+ 3.9518 

+ 1.3457 

0.0 

stable 

AAA S 


2 

+ 1.4070 

-1.1895 

0.0 

s tab 1 e 



I 0 

+ 1.3906 

I -1.0000 

±0.6124 

unstable 


1 

J 

3 

-0.0775 

i -0.1562 

j — 

0.0 

stable 
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Table 7 (Continued) 


Case E fi ' = - 0. 1 3 


Sub-case 


a 

u .i 



Remarks 

Ei 

+3.0 


+ 

1.4877 

— 






B 

+ 

0.0863 

+0.1707 

0.0 

stable 

e 2 

0. 0 

B 

4 * 

0.5269 

+0.6643 

0.0 

stable 



i 

- 

0.9573 

- 

- 


E 3 

-1.5 

H 

+ 

1.6337 

+1.0286 

0.0 

stable 



2 

- 

0. 1397 

-0.6608 

0.0 

uns table 



3 

- 

0.1655 

-0.3679 

0.0 

stable 




- 

0.4923 

- 

- 


E 4 

-3. 0 

1 

+ 

3. 7379 

+ 1.3243 

0. 0 

stable 



y 

+ 

1. 2412 

-1.1618 

0.0 

unstable 



a 

+ 

1.0977 

— 

_ 




K 

- 

0.0806 

-0.1625 

0.0 

stable 

e 5 

-8.0 


+ 

18.5553 

+2. 0463 

0.0 

stable 



2 

+ 

14.5236 

-1. 9847 

0.0 

stable 



0 

+14. 5227 

-1. 9231 

±0. 4865 

uns table 


i 

3 

- 

0. 0302 

-0.0611 

0.0 

s t ab 1 e 
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§8 Discussions 


1. As mentioned before in §6 the expansions of solutions are so 
arranged that they are always convergent. But it is still a problem to 
find out the most efficient expansions for the solutions; in other words, 
to rearrange the series when it is too slowly convergent. This situa- 
tion occurs when q tends to the unity. In such a case we have another 
type of expansion by exchanging the real period and imaginary period 
with each other, so that for the extreme case of q = 1 we have hyper- 
bolic functions instead of circular functions. 

In any way, in order to have suitable expansions of the solutions it 
is necessary to rearrange them into a different form appropriate to 
each case, respectively. 

Also neglected is a proper method to calculate semi-periods 
(a.', oj" i, <uj, ora^'i ) in terms of given constants /3' , yandu. These 
calculations are closely related with the evaluation of q (if A > 0) or 
qj (if A < 0). The numerical processes to find p and X are also 

omitted here. 

These practical problems are, of course, important if we wish to 
obtain the solutions in detail but they are so complicated that their 
discussion will be postponed. 

2. The adoption of non-dimensional analysis in units of the proper 
powers of yS is of course optional. We can treat the calculation without 
any non-dimensional analysis; or we may have another type of non- 
dimensional analysis in terms of fi' or y for example. However, the 
process which we treated in this paper has a slight advantage; firstly, 
in doing so /3' is an absolute constant depending only on the coefficients 
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of harmonics but independent of the initial conditon; secondly, we may treat 
both the cases ft' > 0 or ft' < 0, in the same framework, for example, we can 
draw a diagram such as Fig. 1 on a single sheet; thirdly, the difference 
between ft> 0 and ft < 0 can be neglected automatically by always selecting 
the sign of^or^ Such that , yCft\ > 0. This is caused by the 

fact that the difference between ft > 0 and ft < 0 remains only on the sign 
of hand k but there is no difference on the nature of the critical character, 

etc . 

3. We have assumed J n to be of n-th order of magnitude. But if those 
assumptions are broken down, how shall the situation change? We may 
consider this breakdown as follows within our framework: ft 0, in 
other words a non-effective case of the anti-asymmetrical terms J 3 and J 5 
(/3',-oo). In this case, the range h„ and k, for the critical nature (the 
transition from librational to revolutional) is not restricted in a limited 
area according as y, - - oo , but h and k themselves are limited. There- 
fore we may tend ft - 0 without any essential difficulty in our theory. 

On the other hand if ft' - 0, keeping ft constant the range of h, and 
^ as well as h and k themselves for the critical nature are not restricted 
as 7 , _ _ oo. Therefore, another type of theory would be required which 
could treat not only the small eccentricity case but also the moderate 
eccentricity case. In this case where the anti-asymmetrical terms pre- 
vail over the symmetrical terms, the general feature, which is constructed 
here by adding the anti-asymmetrical terms to the symmetrical terms or 
at most by considering the same order contributions in some meaning from 
both sources, would be broken down. Fortunately, since ft\ for the earth 
is approximately -0.13, it does not seem to necessitate any new theory 
for this respect. 
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4. A higher order theory beyond the one discussed here could be 

developed along the line which we gave in a previous paper 

(1962). However we must keep in mind the restriction made there 

that the leading terms which are, for example, the ones discussed here, 

should prevail to the remaining terms. In other words, if the leading 

terms are so small that the neglected terms play an important role in 

the behavior of the solution, then such an approximation process would 

be broken down. In addition, in the present paper, we have taken into 

account two parameters such as /3 and/3', the former is assumed to be 

third order and the latter is of second order; namely the non-dimensional 

quantity /3 = /3'//3 2/3 is assumed to be far from zero or from the infinity. 

In the previous paper, it is assumed that, for example, J 2 , J 2 J 4 and 

T are all of the third order of J 2 ~ J 4 /J 2 : this means that ’ lf J 2 + J 4 /J 2 
J 6 

tends to zero, J 2 3 etc. should tend to zero as well, in order to be able to 
apply the general treatment described there, otherwise the leading terms 
are so small that the series obtained along the theory would be divergent. 
A similar situation occurs in the present theory in the aforementioned 
sub-section 3. This is also true m the case of the higher order 
theories than discussed. At any rate it is important to consider first the 

most significant parts ox the terms. 

The difference between the previous theory and the present one 
exists in the fact that here the anti-asymmetrical terms are given 
the same importance as of the symmetrical terms by assuming that 
t is of the n -th order of magnitude in the case of small eccentricity. 

J n 
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§9 Conclusion 

In this paper we tried to solve the equations of motion of a 
close satellite near the critical inclination in the case of small eccen- 
tricity under the influence of the potential given by 


V ^ 


i± 

r 



(sin S ) . 


Assumed are as follows: the coefficient of zonal harmonics J n is of 

the n-th order of magnitude and the eccentricity is of the first order 
of magnitude* The terms which are taken in the Hamiltonian are of the 
sixth order of magnitude in this respect; namely, the terms which have 
e 4 J 2 , e 2 J 4 , e 3 J 3 , eJ 2 J 3 , eJ 5 etc. as factor are included* It is also 
noted that the lower order terms do not enter the Hamiltonian except for 
unnecessary constant terms* 

In §6 the explicit solutions are given with classification according 
to the interrelation of fij, y* and u, (as for the definition of these - 
quantities, see Eqs. (4.4), (2.12) and (3.4)), where /3\ strictly depends 
on the coefficients of zonal harmonics J , y* depends on the projected 
angular momentum to the equatorial plane as well, and u„ is the energy 
constant in units of j3 A/3 . This classification is made in order to have 
real expressions for the solutions. Three parameters involved are 
expressed by the non-dime nsional quantities so that one may have the 
actual solutions by factoring a proper power of y3 1/3 , which is assumed 
as of the first order magnitude. 

For the earth /3 1/3 is -t0.0528/a Q ; accordingly, if the eccentricity 
is confined within some range, say 0.05 or 0.10, then the theory devel- 
oped here is applicable. 
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In any way xt xs noted that the transition from the symmetrical term 
prevailing case to the anti- symmetrical prevailing case withxn the 
approximation adopted here occurs at (between Case II and Case III, 
where a great change especially in characteristics occurs)'. 

P[ = P'/p' 1/i = - \- 

This quantity for the earth is approximately -0.13, therefore, 
earth- s potential is in the anti-asymmetrical term prevailing case in 

this respect. 
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2.0 

Fig 2 Characteristics in h* and k*— plane, 
1 l 

h*=h//3 3 , k*=k//3 3 

(3'*=— 3.0,7*=+2.0 (Case i i, or Case 1) 



Fig 3 |3'*=— 3.0, 7*=+1.0 
(Case I 2 . or Case 2') 



Fig 4 (S'*=— 3.0, 7 *=0.0 (Case I 3 . or Case 3') 















Fig 12 0'*=— 0.25. 7*=+2.0 

(Case llli, or Case 1) 




Fig 14 (3'*— —0.25, 7 *— — 1.3 
(Case 1113, or Case 3 ") 










Fig 21 /?'*=— 0.13, 7*=— 8.0 
(Case Es, or Case 5) 



